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By combining classical molecular dynamics simulations and density functional theory total en- 
ergy calculations, we study the possibility of doping graphene with B/N atoms using low-energy 
ion irradiation. Our simulations show that the optimum irradiation energy is 50 eV with substi- 
tution probabilities of 55% for N and 40% for B. We further estimate probabilities for different 
defect configurations to appear under B/N ion irradiation. We analyze the processes responsible 
for defect production and report an effective swift chemical sputtering mechanism for N irradiation 
at low energies (^125 eV) which leads to production of single vacancies. Our results show that ion 
irradiation is a promising method for creating hybrid C-B/N structures for future applications in 
the realm of nanoelectronics. 

PACS numbers: 81.05.ue, 61.72.U-, 83.10.Mj 



I. INTRODUCTION 

Graphene [1] is a two-dimensional all-carbon nanos- 
tructure consisting of sp^-bonded carbon atoms arranged 
in a honeycomb-like lattice. Similar to carbon nan- 
otubes [2], which can be imagined as graphene sheets 
rolled up to form seamless tubes, it has attracted a con- 
siderable amount of attention, especially due to its unique 
electronic properties Q . Pristine graphene is a zero band 
gap material with a high electron mobility, but intro- 
duction of impurities and structural defects [3, 4] can 
lead to band gap opening and Fermi level shifting. This 
makes graphene an interesting candidate for many ap- 
plications in nanoelectronics. However, despite recent 
breakthroughs in graphene production [HI, Q , one of the 
remaining problems on the way towards graphene-based 
electronic devices is the lack of control over the electronic 
properties of the as-grown samples. 

In today's semiconductor industry, a standard way to 
modify the electronic properties of a material is doping 
with foreign atom species. For an all-carbon material, 
such as graphene, the natural dopants are boron and ni- 
trogen, which posses one electron less and more than car- 
bon, respectively. Since they are also roughly similar in 
size to carbon, they can be incorporated into substitu- 
tional positions in carbon structures by replacing exactly 
one carbon atom, or take the vacant site at the edges 
of graphene flakes 04111 • It has been shown [l2| than 
B/N-doped graphene displays p/n-type behavior. Not 
surprisingly, several theoretical studies on hybrid C-B/N 
systems [13-.:18] have recently been carried out. Particu- 
lar attention has been paid to graphene nanoribbons plol - 
[245. Although it is possible to introduce B/N atoms into 
graphene during synthesis (ol-ITTl [25l - [28| , not many stud- 
ies have investigated the possibility of selectively intro- 
ducing B/N impurities into graphene after growth, while 



post-synthesis doping may be an alternative way to cre- 
ate graphene-based functional materials. It was demon- 
strated in a recent study that B can be incorporated 
into graphene structure by selectively exposing different 
sides of a graphene monolayer to different atomic envi- 
ronments [2y]. Unfortunately, it is not clear how this 
could be carried out in practice. While ion irradiation 
is routinely used nowadays when manufacturing conven- 
tional semi-conductor devices, one can expect that it can 
be used in a similar manner to introduce B/N atoms into 
graphene. However, graphene, an atomically thin target, 
has a significantly different response to ion irradiation as 
compared to traditional three-dimensional materials |30| . 
In a recent experimental study [31], N+ ion irradiation 
was used to facilitate N-doping of graphene. Due to a 
high irradiation energy (30 keV), very few if any ions 
were directly incorporated to the graphene sheet. In- 
stead, the nitrogen atoms were later introduced during 
annealing in a NH3 atmosphere. Clearly, the optimiza- 
tion of the implantation process requires full microscopic 
understanding of the response of graphene to ion irradi- 
ation. 

In this study, we employ classical molecular dynamics 
(MD) simulations and density functional theory (DFT) 
total energy calculations to study the efficiency of ion 
beam doping of graphene as a function of ion energy. 
We show that at energies close to 50 eV, it is possible 
to incorporate as much as 55% of N and almost 40% 
of B impinging ions to the substitutional dopant site in 
graphene. The most typical defects which appear un- 
der low energy irradiation include single and double va- 
cancies. In the preferred range of doping energies, the 
probability for introducing additional defects to graphene 
during ion beam doping remains below 37%. We also re- 
port a swift chemical sputtering mechanism (especially 
for N ions) for ion energies below 150 eV, which along 



2 



with the bahistic knock-out, leads to production of sin- 
gle vacancies. This effect may have an important role in 
defect creation during experimental ion beam doping of 
graphene, and possibly other carbon materials. 



II. METHODS 

To study B/N ion implantation into graphene, we fol- 
lowed the same approach as in our earlier work on ion ir- 
radiation of carbon nanotubes (32l-[35j and graphene ^] . 
The interactions between the atoms in the system under 
study were described by the bond order potentials devel- 
oped by Matsunaga et al. [36] for C-B/N and by Bren- 
ner et al. js^ for C-C interactions (without the bond 
conjugation term). Heat dissipation at the boundaries 
of the system was modeled with the Berendsen thermo- 
stat [38] with various time constants. However, for the 
used system size and rather low irradiation energies (up 
to 4 keV) included in this study, the choice of the time 
constant had only a minor effect on the results. Our 
graphene sheet consisted of a 16 x 16 supercell with a 
total of 512 atoms. During the simulations, we directed 
the B/N ion towards graphene in perpendicular direc- 
tion with energies between 10 eV and 4 keV. 1000 im- 
pact points were selected randomly within the irreducible 
area of the graphene lattice to ensure statistically correct 
sampling of the irradiated structure. The total number 
of independent simulations was nearly 1 000 000. 

We repeated the simulations with the Tersoff C-C po- 
tential [39] to get an insight on how sensitive the results 
are with regard to the choice of the potential. No signifi- 
cant differences were observed. Moreover, since the Bren- 
ner potential gives a displacement threshold energy (min- 
imum kinetic energy required to sputter a carbon atom) 
which is very close to that obtained with dynamical DFT 
calculations (22.20 eV as compared to 22.03 eV [40|), we 
are confident that the defect production upon irradia- 
tion is at least qualitatively correctly described within 
our simulation model. The same interaction models were 
also used in our earlier simulations regarding ion implan- 
tation into carbon nanotubes [32, 33]. For these systems, 
our theoretical predictions on nitrogen implantation have 
been later experimentally corroborated [41]. The charge 
of the ion was not explicitly taken into account in our 
simulations, since it is known to have only a minor role 
in the defect production during ion irradiation of carbon 
nanostructures [42|. However, we refer to the projectile 
as an "ion" to facilitate the direct juxtaposition with the 
experimental studies. 

As the C-B/N interactions in our model have been 
tested much less than the C-C interaction, we employed 
DFT total energy calculations to ensure that the result- 
ing defect structures involving B/N atoms, as given by 
the MD simulations, are consistent with ah initio results. 
Our DFT calculations were carried out with the VASP 
simulation package [43|, [44] using projector augmented 
wave potentials [45| to describe core electrons, and the 




Figure 1. (Color online) Typical defect structures in graphene 
observed after low-energy irradiation with B and N ions, as 
calculated with MD method, followed by structure relaxation 
using DFT calculations: (a) perfect substitution of one carbon 
atom with the dopant, (b) dopant- vacancy-complex structure 
in which the dopant atom is bonded to two carbon atoms, and 
(c) dopant atom in the bridge-configuration as an adatom on 
top of a perfect graphene lattice. 

generalized gradient approximation [46| for exchange and 
correlation. Kinetic energy cutoff for the plane waves was 
set to 500 eV, and all structures were relaxed until atomic 
forces were below 0.01 eV/A. Our initial structure con- 
sisted of 200 carbon atoms, so that the finite size effects 
are expected to be small. The Brillouin zone sampling 
scheme of Monkhorst-Pack [4^ with up to 5 x 5 x 1 mesh 
was used to generate the /c-points. 



III. RESULTS AND DISCUSSION 

We started our simulations by calculating the probabil- 
ities for an irradiation event to produce the most typical 
defect structures (some of which are shown in Fig. [1]), 
as functions of ion energy for both B and N ions. The 
most prolific outcomes were substitution (exactly one C 
replaced by the B/N), substitution with a neighboring 
vacancy, single vacancy (SV) and double vacancy (DV). 
At low energies, we also frequently observed the dopant 
atom attached to the perfect graphene structure as an 
adatom. A schematic of the simulation setup as well as 
the defect creation probabilities of these structures are 
presented in Fig. [2] for those defects involving the dopant 
atom (presented in Fig.[T]), in Fig. [3] for SV and in Fig. 2] 
for DV. 



A. Incorporation of the dopant atom into graphene 

The highest observed probability (55%) for any of the 
structures (except for adatom deposition at very low en- 
ergies which naturally gives a probability of 100%) was 
seen for N substitution (Fig.[T]3) at energies close to 50 eV 
(Fig.[2]3). For substitution, it is intuitively clear that the 
probability should display a single peak at these ener- 
gies: although the displacement threshold for a carbon 
atom is about 20 eV in a head-on collision, on average 
the ion must have a somewhat higher energy in order to 
displace the atom, as a nearly head-on collision is un- 
likely. For B substitution the maximum is slightly be- 
low 40% (Fig. and it appears at a slightly lower 
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Figure 2. (Color online) (a) Schematic representation of the 
simulation setup. Probabilities for different defect configu- 
rations as functions of the ion energy: (b) perfect substitu- 
tion, (c) dopant- vacancy-complex structure, (d) dopant as an 
adatom on top of the graphene sheet. 



energy than for N. The lower energy for B ion is a results 
of the difference between the momentum transfer from 
B/N ions to the target atom while the lower probability 
arises from the smaller displacement cross section for the 
B ion. By comparing to noble gas ion irradiation results 
for single vacancy production in graphene ^] (substi- 
tution occurs when the ion is stopped after a success- 
ful sputtering of one target atom), it is easy to see that 
these features are indeed caused by the mass difference 
between the ions rather than chemical effects. For ener- 
gies higher than 0.5 keV, very few ions were incorporated 
into the graphene atomic network. Instead, at high ener- 
gies (where chemical effects are not important, as shown 
below) relatively larger vacancy-type defects will be pro- 
duced [30|. 

For nanotubes [32|, [33|, the highest observed proba- 
bility for B and N substitutions were similar to those 
found for graphene, although the exact values were 
slightly lower /higher for N/B. For both ions, the peaks 
in graphene appear at lower energies and are significantly 
narrower. This can be understood based on the geomet- 
rical differences between nanotubes and graphene. As a 
crude approximation, one can think of the nanotube as 
a parallelepiped with two walls perpendicular to the ion 
beam direction and two walls parallel to it. Therefore, 
it is possible for the ion to be slowed down by the first 
wall before the displacement takes place. This naturally 
leads to substitution events also at energies higher than 
what is observed for graphene, a two-dimensional system. 
What is perhaps more surprising is that the probabilities 
themselves are not much different. A priori one could 
have expected the nanotube structure to yield roughly 
twice the number of substitutions due to the two walls 
parallel to the beam. However, since it also takes en- 



ergy to penetrate the first wall, this is manifested only at 
energies above those for the peaks observed in graphene. 

Another type of observed substitutions is a result of 
an irradiation event in which the ion has first displaced 
two carbon atoms from graphene and then remains in 
the place of one of the displaced atoms in a pyridine- 
like configuration [48], as shown in Figdjb). We never 
observed pyrrolene-like structures [48] in which the ni- 
trogen/boron atom occupies a site in a five-membered 
atomic ring. This situation may be more probable in 
nanotubes [26] where such configurations are assisted by 
the curvature. Since it was not seen in the nanotube 
ion irradiation simulations |32|, [33| either, it may also be 
a feature which only occurs when nitrogen dopants are 
introduced during synthesis. 

Since two atoms must be displaced for a dopant- 
vacancy-complex to be created, it correspondingly oc- 
curs at higher energies than the perfect substitution. The 
probability for this process is presented in Fig. [2l3 as a 
function of the ion energy. Surprisingly, the probabili- 
ties are significantly different for B and N ions, with a 
probability for N being larger by more than two orders 
of magnitude. However, even in this case the probability 
maximum is only ~ 3%. Again, the difference is mostly 
due to the mass difference between the B and N ions. 
For the heavier N, with increasing energies the scatter- 
ing process leads to displacements of the carbon atoms 
in the directions perpendicular to the initial ion direc- 
tion, as described in Ref. [30], while the ion only slightly 
deviates from the original direction. Hence, the defect 
is produced by simultaneous displacement of two carbon 
atoms while the ion itself stops at the graphene sheet. 
However, for B, which is lighter than C, the situation is 
opposite; with higher energies the ion itself is scattered in 
a direction parallel to the plane where it will be scattered 
again by another C atom and thus eventually leaves the 
structure. 

The probability for creating the dopant-vacancy- 
complex in nanotubes [32|, [33] was similar to graphene 
values for N but significantly higher for B (at optimal en- 
ergies). This is due to the fact that the projected atomic 
density in the direction parallel to the initial ion trajec- 
tory is in some areas much higher for nanotubes than 
for graphene due to curvature. Therefore, even though 
the displacement cross section for the B ion with respect 
to target atoms is relatively small at energies required 
to displace two atoms, simultaneous displacement of an 
atom pair is still more likely than in the case of graphene. 

As mentioned above, we checked all observed final 
structures with DFT calculations to ensure their stability 
beyond the description provided by the bond-order po- 
tential. Indeed, both of the substitution cases remained 
stable during the structure optimization (the DFT- 
optimized structures are shown in Fig.[T]). The formation 
energy for the perfect substitution is ca. ~ 4.42 eV lower 
than that for the pyridine-like two-coordinated configu- 
ration. Therefore, one can expect that in the case of 
migrating carbon adatoms on the graphene surface, the 
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pyridine- like configurations tend to get filled with carbon 
atoms to form a structure similar to perfect substitution. 
As has been shown previously [33|, the formation ener- 
gies of the substitutional impurities in graphene are a few 
electron-volts with a somewhat higher formation energy 
for boron substitution. 



B. Creation of single vacancies 

In addition to substitutional impurities, ion- 
implantation should inevitably result in the formation of 
various irradiation-induced defects in graphene [42|, |49| . 
Among these, single vacancies have the highest prob- 
ability to appear. In the case of N irradiation, it 
reaches its maximum at ca. 125 eV with a value of 
55% (Fig. [3ti). This very narrow peak overlaps with a 
broad peak with its maximum (35%) at ca. 400 eV. 
Surprisingly, for B, only a slight effect of the first peak 
is seen as a shoulder of the second peak. In this case 
the peak has a maximum of 35% at 180 eV (Fig. 
The narrower low-energy peak results from a chemical 
interaction between the ion and one of the target carbon 
atoms. From a case-by-case analysis of the distance 
between the ion and the sputtered atom as a function of 
time, it is clear that the attractive interaction between 
the projectile and recoil atom is responsible for the 
production of single vacancies at energies within this 
peak via a mechanism where the ion pulls one carbon 
atom with itself while penetrating the sheet, as sketched 
in Fig. [SJd. This effectively results in the formation of a 
dimer, which may eventually be broken, as also sketched 
in Fig. [3]d. 

In Fig. [3t,d we show the effect of the attraction by 
subtracting those cases where the ion and a carbon atom 
form a dimer during the sputtering process from the over- 
all probability for the single vacancy creation. For B, this 
process completely accounts for the shoulder of the sin- 
gle vacancy probability peak. However, for N, it only ac- 
counts for the lower energy end of the narrow peak. This 
effect can also be directly explained by comparing the 
masses of the ions, and taking into account that both N 
and B have an attractive interaction with carbon atoms. 
Because B is lighter than C, it is unlikely that it will es- 
cape the carbon atom after successfully displacing it at 
the energy range where the attraction plays a role. Thus, 
the dimer remains stable after the sputtering event. How- 
ever, for the heavier N, the ion can effectively pull one 
carbon atom after itself but still retain enough kinetic 
energy to escape the atom thus breaking the dimer. By 
carefully comparing cases with exactly the same impact 
parameters, except for a slight change in kinetic energy 
of the ion, in which the lower energy produced a single 
vacancy and the higher energy did not, we concluded that 
this swift chemical sputtering process (in principle similar 
to what has been observed in Ref. [50] for hydrogen ions) 
indeed completely explains the observed narrow peak in 
the single vacancy production probabilities for both N 
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Figure 3. (Color online) Probabilities for creating single va- 
cancies in graphene under low-energy B and N irradiation (a) . 
Schematic presentation of the swift chemical sputtering pro- 
cess involving chemical interaction between the incoming ion 
and the recoil carbon atom is shown in panel b. In panels 
c and d these probabilities are shown for N and B ions (re- 
spectively) with and without those cases in which a dimer is 
formed between the incoming ion and the sputtered carbon 
atom. 



and B ions. 

These results are in line with Ar irradiation simula- 
tions carried out for nanotubes [34], as well as with the 
noble gas irradiation studies for graphene [30| . If one ig- 
nores the chemical sputtering which does not occur for 
the noble gases, the single vacancy probabilities for B/N 
appear at energies between those for He and Ne in the 
order of increasing mass, as one would expect. Also the 
functional form is similar to that observed for He and Ne 
irradiation, even reproducing the broadening of the prob- 
ability peak with increasing mass of the projectile (3Q| . 
Therefore, also this comparison shows that the additional 
narrow peak must be caused by a chemical effect. 



C. Creation of double vacancies 

The last of the often observed defects is DV. The high- 
est probability (16%) for this defect was found upon N 
irradiation at energies close to 110 eV (Fig. |4^). For B, 
the corresponding peak maximum occurred at 70 eV with 
a value of 6% (Fig.|4)3). In addition to this first peak, the 
probability for B irradiation contains two other, much 
broader, peaks. The first of these (#2) occurs with a 
probability of less than 1% at energies between 0.14 keV 
and 1.4 keV. The second one (#3) appears with a sim- 
ilar probability maximum at energies closer to 2 keV. 
This peak extends to energies higher than those included 
in this study. For N, only the second one of these high 
energy peaks (#3) appears. From the three peaks, each 
of which occur at the characteristic energy range, it's 
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Figure 4. (Color online) Probabilities for creating double va- 
cancies in graphene under impacts of B and N ions. The 
low-energy limit (a) and the probabilities at energies exceed- 
ing 200 eV. Note the difference in the scale of the y axis. The 
numbers in the plot indicate the peaks as discussed in the 
text. 




Figure 5. (Color online) Example trajectories for the different 
double vacancy creation mechanisms, (a) A direct impact by 
the ion with two carbon atoms. The ion strikes in the middle 
of the bond joining these atoms and causes the sputtering of 
both, (b) B ion, which is lighter than C, gets scattered in the 
in-plane direction of graphene so that it has enough kinetic 
energy to cause sputtering of a second carbon atom, (c) Ion 
displaces one carbon atom which then causes another carbon 
atom to sputter. 



easy to predict that there are three mechanisms for the 
double vacancy production. These three mechanisms are 
illustrated in Fig. [5l 

The lowest energy peak is due to a direct impact by 
the ion with two carbon atoms when the ion strikes in the 
middle of the bond between these atoms (Fig. [5^). For 
energies exceeding the value corresponding to the maxi- 
mum in the curve, the scattering cross section for the ion 
and a target atom becomes so small that it is impossi- 
ble for the ion to displace two atoms simultaneously. For 
the heavier N ion, the correspondingly larger cross sec- 
tion extends the probability peak over the energy range 
observed for B. Due to the larger displacement cross sec- 
tion also the peak maximum appears with a higher value 
for N. This also explains the difference observed for the 
ions in the case of vacancy-dopant-complexes, which are 
directly linked to double vacancies. 

The second mechanism (Fig. [SJd), which does not oc- 
cur for N, arises from the mass difference between the 
B ion and carbon atoms. Because B is lighter than C, 
at the intermediate energies it gets scattered in the in- 
plane direction so that it has enough kinetic energy left 
to sputter a second carbon atom. In the case of boron, 
this mechanism leads to the highest probability for the 



vacancy-dopant-complex defect in the case where the ion 
has enough energy to displace both carbon atoms but be- 
comes trapped in the carbon network. For the heavier N, 
this never occurs since it is rather the carbon atom which 
gets scattered in the in-plane direction. Similar behavior 
is also evident in the noble gas irradiation results Q 
where two peaks are observed for the Ne ion irradiation 
(Ne mass is about 1.4 times larger than that of N), similar 
to the N irradiation in the present study. 

The third of the observed double vacancy creation 
mechanism (Fig.[5fc) involves the displacement of one car- 
bon atom in the in-plane direction by the ion, which then 
causes the sputtering of a second atom. In the case of no- 
ble gas ions heavier than Ne, the two peaks (correspond- 
ing to peaks one and three for B irradiation) overlap so 
that the separation of the peaks becomes impossible ^30| . 

D. Dopant as an adatom 

Obviously, for ion energies below the threshold for cre- 
ating a single vacancy, the dopant will be attached to 
the pristine graphene structure as an adatom (if it does 
not go through the middle of the hexagon). Depending 
on the impact point and the energy of the projectile, it 
is also possible that the ion gets bounced back from the 
graphene sheet, which leads to a non- monotonous prob- 
ability curve as a function of the ion energy (Fig. [2]i). 
Although it is rather unlikely that the adatom configura- 
tion would spontaneously transform to a substitutional 
defect, it is possible that the adatoms will take the po- 
sition of a missing atom in a single or double vacancy 
created during the irradiation, since the migration bar- 
riers for B/N adatoms are in the range of 0.1-1.1 eV on 
graphenic surfaces 0, [5l| . 

From a simple formation energy consideration at the 
thermodynamical equilibrium, it would also be beneficial 
for nitrogen dimers to break in order to fill in single va- 
cancies in graphene (according to DFT calculations, this 
would lead to a significant energy gain of more than 6 eV 
per nitrogen atom). This could occur during annealing 
in a N2 atmosphere. However, this process would re- 
quire the breakage of the N dimer which is among the 
strongest molecules in nature. Recent experiments in- 
deed showed [3l| that annealing in a N2 atmosphere does 
not lead to N-doping of defective graphene. However, 
using a molecular environment in which the nitrogen (or 
boron) atoms are less strongly bound, such as NH3, will 
provide dopant atoms to fill the vacancies, as was demon- 
strated in Ref. [31] . 

E. Relative probability for substitution 

Naturally, one of the critical factors in using ion ir- 
radiation to dope a nano-object such as graphene, with 
plenty of free space for the target atoms to escape, is 
the amount of other defects created during the implan- 
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Figure 6. (Color online) Relative probability for B/N sub- 
stitution as compared to creation of any other defects as a 
function of ion energy. 

tation. Luckily, the production of single vacancies, and 
especially more complex defects, occurs mainly at ener- 
gies higher than that of the highest substitution proba- 
bility. The total probability for creating any other defect 
during nitrogen ion irradiation at 40 eV is 29%. For 
boron, the corresponding value at 100 eV is 27%. These 
energies correspond to 10% substitution probabilities at 
the lower and higher end for nitrogen and boron, respec- 
tively. At the most efficient doping energy (50 eV) for 
N, the probability for creating any other defects is 32%. 
For B at this energy, in 36% of the events are other de- 
fects created. These relative probabilities, defined as the 
ratio of sum of probabilities for the different substitution 
cases (S^P/^^®*) to the sum of probabilities for all de- 
fects (S^i^"°«^bst ^ 5]./>subst)^ 3/1^ substitution are 

shown in Fig. [6] as a function of the irradiation energy. 
An obvious way to overcome any problems arising from 
creating a substantial amount of single vacancies during 
the substitution is to combine irradiation with two ener- 
gies; one at the substitution maximum and another one 
below the single vacancy creation threshold. This would 
allow deposition of dopant adatoms which could fill in 
the created vacancies. Although we are unaware of any 
experiments where such a method would have been used, 
this is in many ways similar to that of annealing in a suit- 
able molecular atmosphere providing the dopant atoms, 
as presented in Ref. [31j . 



cient doping energies are, according to our simulations, 
close to 50 eV for both ions. The high-probability peaks 
are wide enough to allow doping also at various ener- 
gies both below and above the maximum. For exam- 
ple, substitution probability for nitrogen ion remains 
over 10% for energies between 40 eV and 105 eV, and 
for boron between 40 eV and 100 eV. Assuming a con- 
trol over ion energies and the possibility to use focused 
ion beams, one can envisage fabrication of spatially- non- 
uniform graphene based materials where N and B doped 
areas can be made with a resolution of better than 5 nm. 
As the defect structures which we observed during our 
MD simulations were confirmed to be at least metastable 
by our DFT calculations, and because our earlier theoret- 
ical predictions [s?, '33] on ion implantation into carbon 
nanotubes have been experimentally corroborated [4l|, 
one can expect that our results for graphene give a good 
indication of the outcome of actual experiments. 

As a conclusion, our computational work combining 
classical MD simulations and DFT total energy calcula- 
tions indicate, similar to what has been shown for carbon 
nanotubes, that ion irradiation is a promising method for 
creating hybrid C-B/N structures for future applications 
in the realm of nanoelectronics. The beneficial irradi- 
ation energy proved out to be 50 eV with substitution 
probabilities of 50% for N and 40% for B. At these en- 
ergies the probability for creating other defects is close 
to 30-50%. By lowering the implantation energy below 
the threshold for creating single vacancies or by anneal- 
ing in a suitable molecular environment (such as NH3), 
the vacancies created while doping can be eliminated via 
dopant atoms to reach a significantly higher substitu- 
tion/defect ratio. We also reported a surprisingly effec- 
tive chemical sputtering method for N irradiation at low 
energies (125 eV) which leads to production of single va- 
cancies. 
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